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AN  ANALYSIS  OF  THE  BREAKUP  OF  SATELLITE 
1974-103A  (COSMOS  699) 


INTRODUCTION 

On  17  April  1975  the  Department  of  Defense  Space  Surveillance  System  (NAVSPA- 
SUR,  operated  by  the  Navy,  and  SPADATS,  operated  by  the  Air  Force)  observed  that 
satellite  1974*103 A (Cosmos  699)  had  disintegrated  into  several  fragments.  According 
to  the  NASA  Goddard  Satellite  Situation  Report,  SPADATS  identified  50  objects  associ- 
ated with  this  breakup.  Orbits  for  38  fragments  were  determined  by  NAVSPASUR  in 
the  course  of  their  analysis  of  the  event. 

The  breakup  of  Cosmos  699  occurred  as  the  NRL  Systems  Research  Branch  of  the 
Space  Systems  Division  was  investigating  the  dynamics  of  ensembles  of  particles  such  as 
would  result  from  a satellite  breakup.  The  investigation  had  shown  that  the  application 
of  statistical  mechanics,  as  done  in  the  continuum  theory  of  stellar  dynamics,  was  a particu- 
larly promising  approach.  The  interest  surrounding  the  Cosmos  699  breakup  has  made  it 
a focal  point  for  the  testing  of  the  newly  developed  statistical  theory  of  satellite  breakups. 

This  report  describes  an  analysis  of  the  Cosmos  699  breakup  base  on  the  application 
of  statistical  mechanics  and  using  NAVSPASUR  elements  for  observational  data.  A guid- 
ing principle  of  the  analysis  is  to  contrast  at  every  opportunity  results  obtained  from  the 
individual  elements  and  those  obtained  from  the  statistical  structure  of  the  cloud  as  a 
whole, 


THE  OBSERVATIONAL  DATA 

The  analysis  is  based  on  orbital  elements  of  individual  fragments  which  were  determined 
at  NAVSPASUR  Headquarters,  Dahlgren,  Virginia.  Thirty-eight  element  sets  were  supplied. 

Of  these,  three  were  rejected  because  of  anomalous  mean  motion  rates,  and  two  were  rejected 
because  the  time  of  epoch  was  considered  to  be  either  too  near  or  too  far  from  time  of  break-up. 
Thus  33  element  s.s  were  retained  for  analysis.  The  elements  are  listed  in  Table  1.  Ele- 
ment set  1 is  for  ne  original  COSMOS  699  payload.  This  element  set  provides  the  refer- 
ence or  parent  orbit  for  the  analysis. 

The  NAVSPASUR  differential  corrections  (DC’s)  which  estimated  the  elements  were 
based  on  direction-cosine  measurements  of  the  NAVSPASUR  “fence.”  The  average  num- 
ber of  observations  per  DC  was  ten  over  a time  span  of  3 days.  Rou  jhly  2/3  of  the  obser- 
vations were  single-station  direction-cosine  measurements,  and  roughly  1/3  involved  two  or 
more  stations,  which  permitted  a determination  of  range  by  triangulation.  The  DC  for  each 
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element  set  used  estimated  the  cl?;ssical  elements  and  the  rate  of  change  of  mean  motion 
which  parameterizes  the  atmospheric  drag  perturbation.  The  rms  value  of  the  residuals  for 
the  DC’s  was  typically  500  to  1000  meters. 

For  the  present  analysis  the  element  sets  were  propagated  by  calculating  two-body 
and  secular  perturbative  effects.  This  introduces  an  inconsistency,  because  the  complete 
Brouwer  theory  was  used  in  the  DC  process  at  NAVSPASUR.  However  the  differences 
(certainly  less  than  1 km)  are  insigniflcant  for  the  purposes  of  this  analysis. 

The  elements  were  found  to  fall  into  two  classes.  The  first  class  consisted  of  28  frag- 
ments which  dispersed  simultaneously  from  the  parent  body  at  107** .9083  (3.0°N,  82.9®W). 
The  second  class  consisted  of  five  fragments  which  apparently  fragmented  at  different  times. 
The  second  class  will  be  discussed  further  in  the  concluding  section. 

The  evolution  of  the  positions  of  the  fragments  is  shown  in  Figs.  1 through  8 during  the 
first  1/2  revolution  after  breakup.  The  differences  are  resolved  along  the  unit  vectors 
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Fig.  1 — PotiUon  of  (ngmenU  nUUw  to  th«  payload  whan  v ' Tq  ■ 0.03  radian  (a)  (t) 
plane  (b)  plane,  Pragmenta  1 through  S correepond  to  element  leta  9,  29,  31,  32,  and 
33  rcapectively  (Table  1). 
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Fig.  2 — Position  of  fragments  relative  to  the  payload  when  v - Vg  « 0.48  radian 
(a)  Jt}  pUne  (u)  plane 
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Fl|.  7 — Position  of  fragments  relstiro  to  the  payload  when  v - V0  ■ 2.91  radians 

(a)  {f?  plane  (b)  plane 


Fig.  8 — Position  of  fragments  relative  to  the  payload  when  v - Vq  » 3.40  radians 

(a)  in  plane  (b)  plane 
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•l  = r/  |r| 

63  = r X v/  |r  X v| 

•2  ” *3  ^ *l 


((  component), 
(f  component), 
(t?  component). 


where  r and  v are  the  position  and  velocity  vectors  of  the  parent.  Figures  la  through  8a 
show  the  projection  of  the  relative  positions  onto  the  ^ 17  plane,  and  Figs,  lb  through  8b 
show  the  projections  onto  the  ^ ^ plane.  The  five  fragments  of  the  second  class  are  indi- 
cated by  the  numbered  filled  circles  in  Fig.  la. 


THEORETICAL  DEVELOPMENT 

The  theoretical  basis  for  this  analysis  is  the  application  of  statistical  mechanics  to  study 
ensembles  of  noninteracting  particles  [1].  There  is  however  an  area  of  the  theory  which 
must  be  expanded  before  we  can  proce^  with  the  analysis.  The  inverse  problem  of  deter- 
mining the  age  and  initial  velocity  distribution  has  been  treated  in  detail  only  for  a two- 
dimensional  ellipsoidal  breakup  from  a circular  parent  orbit.  This  must  be  extended  to 
three  dimensions  in  order  to  proceed.  Thus  this  section  is  devoted  to  the  inverse  problem 
for  a slowly  dispersing  three-dimensional  ellipsoidal  breakup.  There  is  no  reason  for  re- 
striction to  circular  parent  orbits  in  the  theoretical  development,  but  the  eccentricity  of 
Cosmos  699  was  sm^l  enough  that  this  simplification  can  be  used  in  the  application. 

Statistical  mechanics  theory  treats  the  ensemble  of  fragments  as  a continuum  of  par- 
ticles described  by  a phase-space  distribution  function  f(q,  p,  t).  if  the  equations  of  mo- 
tion of  the  individual  fragments  are  determined  by  the  Hamiltonian  }{  (q,  p,  0>  then  f 
satisfies  Liouville’s  equation 

(II 

where  (•,  ■)  denotes  the  Poisson  bracket.  The  initial  condition  satisfied  for  a body  disin- 
tegrating at  coordinate  at  t = 0 is 

Aq.  P.  0)  * 5(*»  - Q*)  G (P).  (2) 

where  5(-)  is  the  Dirac  delta  function  and  G(p)  describes  the  initial  distribution  of  momenta. 

The  equations  of  motion  may  be  linearized  about  the  reference  trajectory  if  the  parti- 
cles disperse  slowly  or,  equivalently,  if  the  concern  is  with  the  initial  phase  of  the  breakup. 
Then  the  solution  of  the  equations  of  motion  becomes 


o-ersa- 


As  shown  in  Ref.  1 the  solution  for  f becomes 


f = S(l/_q  + V_p  - q J G(iv:q  + F.p), 
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and  the  spatial  density 


becomes 


P(q,  0 = /rt».  P , 0 dp 

p - det(V:l)  0[(lV.  - Y.VZ^U.)  q + V.lC'q  J, 


(5) 


(6) 


where  the  notation  U{~t)  ^ U_{t),  etc*»  is  used.  We  may  set  q«  * 0,  because  the- breakup 
originates  on  the  parent  orbit. 


The  momentum  distribution  for  an  ellipsoidal  breakup  may  be  written 

0(p)  - const  e’^"  ’ ‘ '*o^- 


(7) 


The  problem  is  to  find  the  mean  momentum  P0  and  the  symmetric  dispersion  matrix  A 
given  the  position  of  the  fragments  relative  to  &e  parent  body  as  a function  of  time.  To 
do  this,  the  spatial  density  function  is  written  as 


with 


X(f)e^ 

(8a) 

(Idet  VJ)-l, 

(8b) 

u • Qu, 

(8c) 

q -Af'Jpo, 

(8d) 

MT  AM, 

(8e) 

w.  - yLvj^c/.. 

(8f) 

and 

i 

The  initial  mean  momentum  is  obtained  immediately  as 

Po  * Mq. 

where  q is  the  mean  of  the  relative  positions  at  time  t.  To  obtain  A,  the  matrix  J is  de- 
fined by 


entjjJ  - J UiUjP  du. 


(9) 


This  matrix  is  the  point  of  contact  with  the  observations  and  is  obtained  by  calculating 
the  second  moments  about  the  mean  of  the  relative  coordinates  of  the  fira^ents.  The 
theoretical  value  of  the  integral  (9),  assuming  the  ellipsoidal  distribution  (7),  then  relates 
A to  the  observables.  If  (8a)  is  substituted  into  (9),  then  [2]  the  result  is 


J = \ x3/2  K(det  0)-i/2  Q.1 


(10) 
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The  determinant  of  (10)  yields  the  relation 


(det  (det 

\ir3/2  k) 


Therefore 


Q = (det  J)i/5 


Finally,  from  (8e), 


A » (Af*l)T  Q Af-l, 


(11) 


(12) 


(13) 


and  this  is  the  expression  relating  ^4  to  the  observable  J via  equation  (12). 


The  expressions  (3)  through  (IS)  are  valid  for  any  linear  dynamical  system.  The  anal* 
ysis  here  of  Cosmos  699  however  will  be  restricted  to  linearization  about  a circular  Keplerian 
orbit.  In  this  case  explicit  expressions  for  and  K are 


2(1  - c) 


Af-l  - i (2(c  - 1)  4s  - 3nt 


and 


K 


n 

Is  l-3nts  + 8(1  - c)l  r 


where  s = sin  nt,  c - cos  nt,  and  n is  the  mean  motion  of  the  parent  orbit.  Equation  (11) 
also  provides  an  algorithm  for  estimating  the  age  of  the  breakup.  From  equation  (8e) 

(det  Q)*i  « det  Af*^. 

Therefore  firom  equation  (11) 

(’•312  ir\3l6 

(detAf-1).  (16) 

The  left-hand  side  of  (16)  is  obtained  directly  from  the  observations  as  a function  of  clock 
time.  The  right-hand  side  is  a known  function  of  time  since  breakup.  The  time  of  breakup 
IS  determined  by  translating  the  left-hand  side  in  time  to  match  the  right-hand  side.  The 
translation  determines  time  since  breakup  and  therefore  determines  the  breakup  time  itself. 
Since  the  reference  orbit  is  known,  knowledge  of  breakup  time  immediately  provides  the 
breakup  position. 
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ANALYSIS  OF  FIRST  MOMENTS 

In  this  section,  which  begins  the  analysis  of  the  Cosmos  699  data,  the  first  moments, 
or  means,  of  the  relative  positions  of  the  fragments  are  analyzed.  According  to  the  theory 
in  its  linear  approximation,  the  center  of  mass  of  the  fragment  ensemble  will  execute  an 
epicycle  in  the  reference  plane.  The  out-of-plane  motion  is  harmonic  and  is  decoupled 
from  the  in-plane  components.  The  observed  epicyclic/harmonic  motion  will  be  us^  to 
determine  both  the  mean  of  the  initial  velocity  components  and  breakup  age.  The  virtue 
of  analyzing  the  fragment  ensemble  as  a whole  as  opposed  to  straightforward  statistics 
derived  from  the  NAVSPASUR  elements  is  that  the  former  minimizes  the  effect  of  error 
that  is  inevitably  present  in  the  orbit-determination  process.  The  velocity  differences  at 
breakup  are  comparable  to  the  error  expected  in  individual  elements.  However  the  extent 
of  the  ensemble  cloud  approximately  IjA  revolution  after  breakup  far  exceeds  the  error 
expected  in  the  elements.  Thus  one  can  obtain  reliable  direct  estimates  of  the  structure 
of  the  fragment  ensemble  as  opposed  to  unreliable  direct  estimates  of  the  velocity  differ- 
ences at  breakup  time.  The  statistical  theory  enables  one  to  translate  the  former  into  infor- 
mation about  conditions  at  breakup  time. 

The  mean  of  the  f component  versus  time  is  shown  in  Fig.  9.  The  data  fit  the  sine 
curve 

f * 14.4  sin(i;  - t?o  - 0.03). 

The  theoretical  expression  for  this  mean  is 

f = (Fo/n)  sin(n(f  - Iq)!. 

The  clock  times  associated  with  the  difference  in  true  anomaly  v - Vq  and  the  mean  motion 
of  the  parent  (n  - 0.001123  rad/s)  yields 

^ = 16.2  m/s 
and 

to  = 107<<.9083. 

The  mean  of  $ versus  the  mean  of  17  is  shown  in  Fig.  10.  The  theoretical  expression 
for  this  curve  is  given  parametrically  by 

? = (^o/n)  sin(n(t  - fg)]  + (2i)oln)  {l  - cosln(t  - tg)!} 
and 

ij  = (-2fg/n)  {1  - costn((  - tg)]}  + (^g/n)  {4  sinlnlf  - tg)]  - 3n(t  - tg)}. 

The  data  are  well  fit  by  the  theoretical  curve  with 


,-r  »<  ''^i'f''^T''^<^i^i'.''^''^’  'V 
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Fig.  10  — Average  ( coordinate  veraua  the  average  t)  coordinate 


Tfo  = 16.6  m/s. 

The  amplitude  of  the  sine  curve  and  the  epicycle  are  quite  sensitive  to  fo  and  7iq  r^pectively. 
Consequently  the  uncertainties  in  these  values  are  less  than  0.2  m/s.  The  value  of  S'q  i^more 
difficult  to  estimate.  The  timing  along  the  epicycle  is  the  parameter  most  sensitive  to  ^'q. 

The  value  given  is  probably  accurate  to  1 m/s. 
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In  Fig.  11  are  shown  the  resolutions  of  the  plane  components  of  velocity  residuals. 
Also  shown  are  the  estimated  values  of  Jq  and  r\Q.  They  are  seen  to  be  consistant.  How- 
ever a straightforward  calculation  of  and‘^0  from  these  data  produces  an  epicycle  which 
fits  the  data  poorly.  Hence  this  is  the  first  example  of  the  advantage  of  the  statistical 
theory  with  regard  to  obervational  error. 


Fig.  11  *-  Projection  of  the  velocity  difference  at  breakup 
time  onto  the  $i)  plane 


ANALYSIS  OF  SECOND  MOMENTS 

The  second  moments  contain  a wealth  of  information.  The  breakup  time  may  be 
infened  from  the  run  of  either  the  dimensions  or  the  orientation  of  fragment  ensemble. 

The  velocity  dispersions  may  be  obtained  by  assuming,  as  here,  an  ellipsoidal  distribution. 
These  dispersions  give  some  indication,  based  solely  on  dynamical  considerations,  of  the 
mechanism  responsible  for  the  breakup.  To  sec  why  this  is  plausible,  one  can  consider 
the  two  extreme  possibilities  of  intentional,  on-board  destruction  and  destruction  from 
the  impact  of  a projectile.  The  former  could  likely  yield  a nearly  isotropic  velocity  dis- 
tribution, and  the  latter  should  yield  a highly  directional  velocity  distribution-  indicative 
of  the  angle  of  approach  of  the  projectile. 

The  velocity  dispersions  will  be  considered  first.  Figure  12  shows  the  reciprocal  square 
roots  \ of  the  eigenvalues  of  matrix  A.  This  matrix  was  calculated  from  equations  (12) 
and  (13),  the  second  moments  of  the  matrix  J having  been  calculated  from  the  relative 
positions  according  to  the  NAVSPASUR  elements.  These  eigenvalues  are  constants  in  the 
theory.  The  variability  evidenced  in  Fig.  12  may  be  attributed  to  observational  error  and 
the  departure  of  the  actual  spatial  density  from  the  assumed  ellipsoidal  form.  The  larg- 
est of  the  X,-  varies  the  least  and  may  be  estimated  to  be  4.25  with  an  uncertainty  of  ±0.1 
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Fig.  12  — Reciprocal  square  roots  of  the  eigenvalue  of  the  calculated 
velocity  dispersion  matrix  versus  the  true  anomaly  difference 


(relative  uncertainty  2.3%).  The  other  two  may  be  estimated  to  be  1.7  ± 0.3  and  1.20  ± 

0.2  (relative  uncertainty  approximately  17%  in  both).  This  establishes  the  shape  of  the 
velocity  ellipsoial  to  be  roughly  a prolate  spheroid  with  aspect  ratio  2.8:1.  The  orienta- 
tion of  the  spheroid  is  determined  by  the  eigenvector  associated  with  the  largest  X,-.  Fig- 
ure 13  shows  the  orientation  angles  of  thatjigenvector  in  die  ei62e3  frame.  There  is 
considerable  scatter  about  the  mean  values  6 = 31.1®  and  0 ® 76.5®.  The  scatter  places 
the  direction  of  this  eigenvector  inside  a cone  whose  axis  is  in  the  mean  direction  and  with 
an  apex  angle  of  about  10®.  The  scatter  in  the  other  two  eigenvectors  is  too  large  to  per- 
mit a meaningful  estimate  of  their  directions.  This  may  be  attributed  to  observational 
error  and  to  the  computational  ambiguity  of  finding  the  principal  axes  of  an  ellipsoid  which 
is  nearly  a spheroid.  The  dot  product  of  the  unit  vector  associated  with  the  mean  direction 
of  the  principal  axis  (indicated  by  a cr-'ss  in  Fig.  13),  namely  e = (0.1927,  0.8233,  0.6168), 
and  the  unit  vector  associated  with  the  mean  velocity  increment,  namely  e = (0.1701,  0.7059, 
0.6876),  shows  the  angle  between  the  two  to  be  12.4®.  Thus  the  velocity  ellipsoid  is  ori- 
ented roughly  in  the  direction  of  the  mean  velocity  increment.  This  shows  that  the  Cosmos 
699  breakup  was  highly  directional.  The  mean  velocity  increment  and  the  velocity  ellipsoid 
are  illustrated  in  Fig.  14. 

The  time  of  breakup  follows  from  the  second  moments  from  both  equation  (16)  and 
the  discussion  of  it.  A comparison  of  observed  and  calculated  values  for  equation  (16)  are 
shown  in  Fig.  15.  The  amplitude  is  unimportant  as  far  as  the  breakup  time  is  concerned. 

This  figure  shows  that  the  breakup  time  derived  from  the  first  moments  is  consistent  with 
that  obtained  from  the  run  of  the  dimension  of  the  ensemble.  No  simple  shift  of  origin 
would  yield  a better  fit  than  shown. 


CONCLUSIONS 


First,  some  conclusions  can  be  made  with  regard  to  the  five  outlying  particles  in  Fig. 
la.  The  particle  labeled  5 is  merely  a high-velocity  particle  (Av  170  m/s)  which  does 
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Fig.  16  — Function  which  detenninei  the  time  of  breakup  from 
the  matrix  of  lecond  momenta 


originate  with  the  main  ensemble  but  evolves  much  more  quickly.  The  remaining  four  parti* 

cles  however  appear  to  have  fragmented  at  ctitterent  times.  The  relative  trajectories  of  these 

particles  (projected  onto  the  it]  plane)  are  shown  Fig.  16.  The  positions  of  the  particles  at  | 

breakup  time  are  shown  by  large  fillc-d  circles.  When  particle  1 is  followed  backward  in  * 

time,  it  passes  near  the  parent  exactly  two  revolutions  before  breakup.  If  particles  2,  3, 

and  4 are  followed  forward,  it  is  found  that  2 and  4 pass  near  the  parent  exactly  one  revo* 

lution  after  breakup  and  particle  3 does  so  exactly  3 revolutions  after  breakup.  It  appears 

that  these  particles  fragmented  one  and  three  revolutions  respectively  after  the  main  break* 

up  and  did  not  even  exist  at  breakup  time.  Their  positions  shown  in  Fig.  la  are  pmely 

fictitious.  Another  interesting  observation  is  that  all  four  epicycles  leave  roughly  the  same 

dimensions  (as  do  their  excursions  out  of  reference  plane)  and  that  these  in  turn  are  typical 

of  particles  in  the  main  ensemble.  There  is  little  doubt  that  they  share  a common  fragmen* 

tation  mechanism  with  the  main  ensemble.  The  scenario  of  a breakup  is  summarized  on  a i 

ground  trace  in  Fig.  17. 

Second,  some  comments  can  be  made  with  regard  to  the  effectiveness  of  the  theory. 

In  spite  of  the  facts  that  the  eccentricity  of  the  reference  orbit  is  not  zero  (albeit  small), 
that  drag  was  included  in  the  orbit  calculations  but  not  the  statistical  theory,  and  that  a 
rather  small  sample  size  was  available,  the  theory  works  remarkably  well.  An  examination 
of  Figs,  la  through  8a  quickly  convinces  one  that  the  actual  spati^  density  is  not  ellipsoi* 
dal  either,  because  of  the  apparent  asymmetry  along  the  major  axis.  Thus  higher  order  mo- 
ments would  be  necessary  to  describe  the  original  velocity  distribution  in  more  detail.  How* 
ever,  the  computational  difficulties  encountered  with  the  second  moments  indicates  that  one  f 

must  anticipate  trouble  if  such  an  attempt  is  made.  i 
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